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We study first- and second-order theories of relativistic diffusion coupled to hydrodynamics under 
the approximation, valid at mid-rapidity in the RH1C and LHC, that conserved number densities are 
much smaller than the entropy density. We identify experimentally accessible quantities of interest, 
and show that the first- and second-order theories may lead to radically different evolutions of 
these quantities. In the first-order theory the memory of the initial state is almost completely 
washed out, whereas in the second order theory it is possible that freezeout occurs at a time when 
■ transient dynamics is still on, and the memory of the initial state remains. There are observational 

' consequences which we touch upon. In the first-order theory, and for initial conditions when the 

second-order theory mimics the first-order, one may be able to put a bound on the diffusion constant. 

<n ; 

' Event-to-event fluctuations of conserved quantum numbers in heavy-ion collisions have been espoused as signals of 
the underlying thermodynamics Quantities which have been considered interesting include the baryon number, 
B, the electric charge, Q, and even the strangeness, S, which is conserved in strong interactions. The proton number, 
N p , has been suggested as a proxy for B (2j. Of interest are the net conserved quantum number and its distribution 
over the ensemble of events at the collider. 

In every observation of a collision event, one can find the density of B, Q or S as a function of rapidity. In the 
usual measurement of fluctuations one usually sums such a density over all rapidity within the observational window: 
Qh this corresponds to taking the Fourier coefficient k = of the density profile. However, other Fourier coefficients can 
be easily constructed. In this paper we point out that hydrodynamic evolution of the Fourier coefficients of conserved 
' number densities are interesting in their own right, since they may contain interesting signals not only of the initial 
, conditions but also of hydrodynamic evolution. 

To this end we investigate the coupled evolution of the number densities and the usual hydrodynamic quantities, 
| i.e., the stress tensor, T^, expressed as usual in terms of the field of flow velocity, and the energy density, e, or the 
entropy density, s, and the dissipative parts. In heavy- ion collisions at the RHIC net conserved number densities are 
£■ — ' small compared to the entropy density. This is usually expressed as^the observation that the ratio of the appropriate 
t^- . chemical potential, fx, and the temperature, T, i.e., fi/T, is small p|. This ratio is expected to become even smaller 
at the forthcoming LHC experiment. As a result, the finite /z corrections to the pressure, p, the speed of sound, c s , 
and to e and s are expected to be small. The combination of these quantities along with the number densities obey a 
t-H ■ set of coupled equations. The smallness of \xfT implies that the equations can be linearized in \ijT , and an accurate 
' picture of the hydrodynamics is obtained by considering the diffusion equation coupled to the hydrodynamic flow 
while neglecting the effect of the number densities on the flow. We note that this approximation is made here for 
' convenience. It allows us to understand the essential physics of the situation without much of a sacrifice in accuracy. 
. If higher accuracy is needed, one can easily investigate the fully coupled formalism [J] . 

Previous work had investigated heavy quark diffusion in the heavy-ion induced fireball [f|. The energy loss of the 
heavy-quark was estimated using weak-coupling theory 0], the Boltzmann equation Q, AdS/CFT techniques [H, and 
using Langevin dynamics 0]. A diffusive theory has also been invoked to describe stopping and early time entropy 
production [1 01 ] - Our motivation, and therefore, the theory we develop, is quite different. 

In this paper we consider only longitudinal background fluid flow, i.e., flow in which the z-axis is defined by the 
direction of the incoming nuclei, and the dependence of all quantities on the orthogonal x and y coordinates is 
neglected. This first approximation is expected to be valid in nearly head-on (small impact parameter) collisions of 
nuclei for times shorter than the sound travel time across the fireball, i.e., for r < R/c s , where R is the radius of the 
colliding nuclei. We work with curvilinear coordinates, i.e., the space-time rapidity, r\ = (1/2) log(i + z)/{t — z), and 
the epoch, r = y/t 2 — z 2 . The covariant derivatives, then include Christoffel symbols. We define the local timelike 
derivative, D = u^d^. This allows us to define a vector = (Du^)/S which is a spacelike unit vector orthogonal to 
Mp, where S 2 = —(Du lJ )(Du^). In terms of this new vector we define the spacclikc derivative D = v^d^. The spacelike 
projector is A M „ = — u^Uy. A scalar which appears often is O = d^u^. For longitudinal flow we can parametrize 
the fluid velocity as u T = coshy and Tu n = sinhy with the other two components vanishing. Boost invariant flow 
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corresponds to taking y — everywhere. This is not the same as the fluid being at rest; the choice of the particular 
curvilinear coordinates that we use makes this simple parametrization possible. Then one has = 1/r and S = 
(note that i> M remains finite). The derivative D = d/dr and tD = djdr\. fn all this we follow the notation of [ll|, to 
which we refer readers for further details. 

We examine the kinematics of the number current, n^. One finds the number density as usual, n = n^u^. Then 
one can write = nu^ + w^, where vv^ is the dissipative part of the number current. Clearly, one has v = —v^n^. 
When there is more than one conserved particle number, one can write an analogous decomposition for each conserved 
quantity, with different n and v for each. For each conserved number, the equation of continuity is 

= d^ri 1 = Dn + n<d + Dv + Sv. (1) 

We will assume in this paper that the background fluid undergoes boost invariant longitudinal flow. However, as a 
result of the assumption that fi/T <C 1, there is no contradiction in assuming that is not boost invariant while the 
background flow is. In fact, the physics of diffusion, which is the central object of interest here, would be trivially 
absent if n and v were both independent of r\. 

In ideal fluids there is no dissipative part to any hydrodynamic quantity, so v = 0. The hydrodynamics then lies 
entirely in the continuity equation and its self-similar solution 

Mr,v) = _nM sQthat rpv) = ^ (2) 
dr t n(T ,i]) t 

Any initial number density profile is then propagated in time unchanged in shape but attenuated as 1/r due to the 
geometry of longitudinal flow. We call this phenomenon Bjorken attenuation. 

Following we define the Fourier transform and the power spectrum of the density, respectively, as 

n(r, fc) — / dr)n(T, rf) exp[— ikr)), and P(t, k) — \n(r, k)\ 2 . (3) 



/2vr 

Note that the measure used here is dr], whereas the invariant volume measure contains rdrj. Hence, the Fourier 
coefficient n(r, k = 0) is not the conserved charge, but rn(r, k = 0) is conserved. For an ideal fluid all Fourier 
coefficients of the number density profile also undergo Bjorken attenuation, which, as we see, is a purely geometric 
phenomenon for boost invariant flows. 

For non-ideal fluids the first order constitutive equation for diffusion is Fick's Law. In the local rest frame of the 
fluid this takes the form 

= n* + Vd'n, (4) 

where T> is the diffusion coefficient. In a general frame one can recast this in the form 

= A^n„ - VA^d„n = vv» -V [d^n - u^Dn] , (5) 

using the projector orthogonal to u. When the last form on the right is contracted with u M , one finds that it vanishes 
trivially. A contraction with i; M relates the dissipative part with the spatial derivative of the number density — 

= -v - VDn. (6) 

Substituting this form of Fick's law into the continuity equation and thereby eliminating one finds the relativistic 
version of the diffusion equation in a fluid undergoing longitudinal expansion — 

= Dn - DVDn - VSDn + n6. (7) 

Note that the last two terms are directly related to the hydrodynamic flow. Since the equation is linear, it can be 
solved by Fourier transforming in the spatial variable and considering the evolution of each Fourier mode separately. 
For later convenience we give the name Fick diffusion to phenomena that arise from this equation. 

The relative importance of flow and continuity versus diffusion can be quantified through the dimensionless quantity 
which compares the 2nd and the last terms in eq. ([7]) — 

A 2 

W = -= = ^ sinh 2 A?], (8) 

where A = r sinh At] is an intrinsic length scale in the density profile, corresponding to a scale A77 in rapidity. When 
W <C 1 then diffusive effects dominate over flow; when it is much greater than unity, flow dominates. Which behaviour 
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dominates depends also on how V changes with time. In the high temperature phase of the plasma, where there is 
essentially only one scale, T, one has D a 1/T (in an AdS/CFT computation one typically obtains T>T = l/2n, 
although much smaller and larger results can also be obtained by tuning parameters [12|). However, at sufficiently 
small T (close to T Cl for example) it is possible that T> is controlled by a hadronic length scale, and may change only 
marginally with T . In the high temperature phase, if the background flow is boost invariant and longitudinal, then 
T oc 1/t( 1+c ^/ 4 . In this case W increases with r, leading to a decreased importance of the diffusive term. At lower 
temperature as well, since T> varies little, W increases and the flow terms become more important. If we are concerned 
with a Ar] range such that diffusion dominates at initial times, then after a time Tfi = 2?/sinh 2 An, flow dominates 
over diffusion. If at times less than Tfi diffusion manages to destroy structures in the number density profile, then 
the memory of the initial state can be lost. 

This can be seen in a simple model of a fluid where all transport coefficients are constant and independent of 
parameters such as [i and T . Then, for any longitudinal background flow of such a fluid the diffusion equation (|7|) 
becomes 

dn n T> d 2 n 
dr t t 2 dn 2 ' 

after linearizing in /i/T '. The Fourier coefficients of n obey the equation 



dn(r, k) n(r, k) T>k 2 , , . . . 

0= a ' +^^ + ^n( T ,k). (10) 



The solution can be obtained by quadrature- 



n(r, k) — ti(to, k) ^— ) exp 



Vk 2 f 1 _TQ 



T \ T 



(11) 



The constant mode (k = 0) exhibits Bjorken attenuation, as all Fourier components eventually do, consistent with 
the analysis of W. Behaviour typical of first-order diffusion is the relative suppression of modes with larger k, i.e., 
the monotonic smoothing of any initial density profile. 

It is possible for experiments to give an upper bound for T>. If one sees no structure in the number density profile 
for rapidity separations up to An then one can conclude either that structures on such scales were not present in the 
initial state, or that they were present and were wiped out by diffusion. The time available for diffusion to act is 
bounded by the freezeout epoch, r/. Hence, one gets the limit T> < Tf sinh 2 An. One should be careful about this 
upper limit for T>, if An is very small, i.e., if one observes fluctuations in the number density profile at all scales. 
There is a sum rule between T> and the relaxation time for diffusive processes, tr [H, [l4|: V — c 2 tr. When freezout 
occurs away from a critical point, a vanishing T> implies a vanishing tr, whereas we expect that tr is not zero. Hence, 
if we see fluctuations of number densities at all scales A77, then it is likely that the first order theory breaks down. To 
find what could happen then, we next investigate the second order theory. 

At second order, Fick's law, eq. H]), can be replaced by the form (T3, [bj| 

= (l + T R d t )n l +Vd i n, (12) 

where tr is a relaxation time for number changing processes. This equation is written in the rest frame of the fluid. 
One uses the projector orthogonal to u to write this covariantly. The covariant equation is 

Q = vv» + trA^Duu - V{d^n - u^Dn). (13) 

Projecting parallel to u gives identically. Projecting parallel to v and using the definition Du^ = Sv^ and the 
identity v^Dv^ = gives 

v + t r Dv = -VDn - Sr R n. (14) 

This replaces the constraint (eq. [H]) for v obtained in the first order formalism. Kelly makes the identification 
T>/t r = c 2 , and we retain this intuition in this work. The quantity t r is currently unknown; it could be of the order 
of typical QCD scales, i.e., 1 fm, or it could be significantly smaller in AdS/QCD scenarios, for example, 0.1 fm. 
In the background of boost-invariant longitudinal flows, eqs. |T]) and (fT4"|) reduce to 

dn n 1 dv , dv T> dn , . 

= — + - + - — , and = t r — +u+— — . (15) 
ot t t on ot t or) 
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This is the form that Kelly's second order diffusion equation [13[ takes in a boost-invariant geometry. In general one 
expects the diffusion constant T> and the relaxation time tr to depend on the temperature, T. Since the time-evolution 
of T is obtained by solving the remaining hydrodynamic equations, one gets an explicit time dependence of V and 
tr, and hence one can solve the equations by Fourier transformation. The evolution equations for the Fourier modes 
can be written in the form 



d_ f n(T,k)\ _ _ f 1/r ik/r\ ( n(r,k) 
dr \ Kt, k) J I ic 2 s k/r 1/tr ) { v{t, k) 



(16) 



Since c 2 varies within a bounded region, c 2 s /t falls with r, and at late times v decays as exp(— t/tr), and n decays as 
1/r, thus reducing to Bjorkcn attenuation at late times. The k = mode exhibits Bjorken attenuation at all times. 
For more detailed analysis it is useful to change variables to i) — log(r /tr). The equations then become 

3 (n{d,k)\ _ .,(n(0,k)\ „_{ 1 ik 



W\v&,k))-- M {vV,k))> M -[ic 2 s k e*)' (17) 
Note three regions of Fourier modes — 

1. When k is sufficiently smaller than unity, the off-diagonal terms may be neglected to a good approximation, and 
the problem decouples. The equation for ra(i?, k) is then similar to the ideal case. We do not discuss the k — 
mode in the following. 

2. When k is sufficiently larger than exp(i?), the problem simplifies again. However, the hydrodynamic equations 
are an approximation to microscopic physics, and are valid for long wavelength spatial fluctuations and for 
long-time phenomena. The time scales are expi? > 1, and the limit of k — > oo must be discarded, although the 
equations simplify. 

3. In the remaining region the equations are fully coupled and a more detailed analysis is called for. This follows. 

Since M depends explicitly on a complete solution to the differential equations cannot be obtained by just 
diagonalizing it. First note that when k ^ 0, then for either § ^ or c 2 s ^ 1, one has [M, M'] ^ 0, i.e., M is 
non-normal. As a result, M is not diagonalized by a unitary transformation, and the eigenvectors are not orthogonal. 



For non-normal matrices the concept of pseudospectra can yield powerful results [15j. In this paper, however, we 
present a standard spectral analysis, since very detailed information can be obtained by this route. 
Since the trace and determinant of M are real, the eigenvalues are either both real or both complex, 

A± = \{l + e*) ± y^~l) 2 -Ac 2 k 2 . (18) 

For a given the eigenvalues are real for |fc| < (e 15 — l)/(2c s ). When they are real, then they are both positive since 
the first term of the expression in eq. (fT5)) is greater than e^/2 and the second term is less than e 9 /2. When both 
are complex, then the real parts are positive, since the trace is so. This implies that \n\ 2 and \v\ 2 both must decay at 
long times. However, there may be transient growth. 

One way to show that transient growth may occur in general is to write down the evolution equation for the power 
spectrum 

P(0,fc) = \n($,k)\ 2 = xUx where x=("j and A=(^ °). (19) 



The evolution equation is 



dP(-8,k) dxt dx t / 2 ik\ . n 

- Ax + x 1 A— = -x 1 Mx, where M = ., , (20) 



dtd M dtd ' \-ik 

independent of c 2 and The quantity — x^.Mx/|a;| 2 is called the numerical range of the matrix — M., and is bounded 
by the eigenvalues of —M,, since the matrix is Hermitean. Since the determinant is negative, one of the eigenvalues is 
positive and the other negative. As a result, the numerical range has indefinite sign, and amplification of the power 
spectrum is possible. The vectors x may be parametrized as x^ = x(sin 9e~ 11 ^ , cos 9), where 8 and <j) range over the 
sphere. Then one has 

- x f Mx = -2\x\ 2 sin 9 (sin0 + k cos 9 sin <j>) . (21) 



FIG. 1: The parts of the sphere in which the power spectrum is amplified are indicated by the green patches. They are 
symmetrically placed at the poles and are bounded in the region < 9 < tan -1 k and symmetrically in 7r + tan~ 1 (— k) < 8 < ir. 

At the equator, 9 — ir/2, this quantity is negative. There are zeroes at the poles, 9 = and ir. Other zeroes occur 
along the curves tan 9 = —ksmtfi. This equation describes two closed curves on the sphere each of which passes 
through one of the poles (see Figure [T]). In the area of the sphere inside these curves (not containing the equator) the 
function is positive. Amplification of -P($, k) can take place whenever x passes through the non-empty region where 
x^Mx < 0. 

Since initial conditions change from event to event, the probability of transient growth is measured by the fractional 
area of the sphere where transient growth may occur. The larger the value of k, the closer does the zero curve come 
to the equator, and hence the larger the probability of transient growth becomes. However, large k corresponds 
to smaller Ary, and hence to short range structures in rapidity. Since such structures are efficiently erased in Fick 
diffusion, their observation in a significant fraction of events would serve as a signal of Kelly diffusion. 




FIG. 2: The evolution of the power spectrum, P($, k), and the number profile, rn(i?,ry), for initial conditions containing only a 
discrete set of non- vanishing modes k — 1/8 (black), 1/4 (red), 1/2 (green), 1 (blue), 2 (dark blue) and 4 (gray), which are taken 
with a flat power spectrum at t = 0. We have chosen (? s = 1/3 for this example. The initial conditions on the corresponding 
^(0, k) are assigned at random. The evolution of the power spectrum is shown on the left and of the profile on the right. In 
contrast to the intuition from the Fick theory, the profile may steepen at intermediate times. 

In Figure [5] we show an example of the evolution of number density profiles using eq. (jTTJ) , starting from random 
initial conditions. In terms of Fourier modes it is clear that transient growth may occur for any k. The typical time 
scale of this growth is 1? « 1, i.e., r « erji. However, even at i? w 2, i.e., r w e 2 Tji the evolution does not begin to 
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FIG. 3: Same as Figure [2] but using the first order diffusion equation. The evolution of the power spectrum is governed by eq. 

CD- 

resemble Bjorken attenuation. If tr is in the range of 1/2 to 2 fm, then the evolution of number densities may be 
dominated by transients. Such long time scales are good news, because they allow one to extract information on the 
early stages of the fireball. 

Figure [3] shows the evolution in the Fick theory of P(i9, k) and Tn{d,k) starting from the same initial conditions 
as above. The difference between first order (Fick) and second order (Kelly) theories of diffusion seems to be fairly 
fundamental. The main intuitive understanding of diffusion from the Fick theory is the following: if initial conditions 
set up sharp density gradients, then diffusion always smooths these out on all scales, monotonically, i.e., as time 
evolves the gradients get smoother and smoother. The Kelly theory can violate this intuition. Depending on initial 
conditions, gradients may be transiently amplified, at different scales at different times. The asymptotic final state is 
the same in both theories. However, since freezeout may not occur at asymptotically late times, the number density 
profiles at freezeout could be different. Some feel of the differences between the theories can be obtained by comparing 
Figures O and [3J Note, in particular, that the solution of Kelly's equation contains structure in the number density 
profile at scales of A77 ~ 1 even at times r ~ 7.5tr. In contrast, for the solution of Fick's equation, the scales of A77 
on which structures are seen are much larger. 

N(i?,A77=l) 

8r 
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FIG. 4: The total conserved charge, N, within a window A77 = 1 (between 77 = and r) = 1) for the initial conditions used in 
Figures [2] and [3] The upper line is for the Kelly theory and the lower one for the Fick theory. 

In FigureS] we show the conserved charge, N(§, Arj = 1), in a rapidity window Arj = 1 (between 77 = and 77 = 1) 
as a function of $. In both theories there is a tendency for the charge to become independent of d at late times. This 
happens because diffusion must eventually wipe out all spatial structure in n(t9, 77), and once the profile has flattened, 
the number density is locally conserved. The differences between the two theories are due precisely to the transients 
that we have studied. In the Fick theory transient lifetimes are extremely small even for k < 1, hence the asymptotic 
behaviour sets in very early. This can also be seen from the profiles in Figure [3J In the Kelly theory, on the other 
hand, transients could be long lasting (depending on initial conditions). This is seen in the evolution of the density 
profiles in Figure as well as of iV(t9, A77 = 1) (Figure 2]). 
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T / |n(r / ,k=l/2)|/o- T / |n(T / ,k=l)|/(r r / |n(T / ,k=2)|/cr 



FIG. 5: The distributions of the square root of the power spectrum of eq. (I22f) in the Fick (blue) and the Kelly (red) theories of 
diffusion for various k at 77 corresponding to 1? = 2. The results are from Monte Carlo simulations with 10000 points for each 
simulation. Initial conditions were chosen from a Gaussian distribution with zero mean and unit variance. All other parameters 
correspond to the values used to obtain the results in Figures [2] and [3] The frequency histograms are normalized to have unit 
area. 



Our discussion till now has been geared towards identifying the fundamental differences between the Fick and Kelly 
theories of diffusion. However, the observables that we have discussed may not be the best suited for heavy-ion 
physics. This is because we have looked at the consequences of starting from a given initial condition. The initial 
conditions for the fireball change from event to event, and averages over events are unlikely to show overdense or 
undcrdense regions. At the LHC particle multiplicities are expected to be high enough that one could build up the 
profile function n(rj, rf) at the freezeout epoch, rj, in a single event. Then by observation of a relatively small number 
of events, one could distinguish between Kelly and Fick diffusion. However, other techniques are required when the 
particle multiplicities are smaller. 

It would be simplest to construct observables that depend on the power spectrum at freezeout. In an event with 
Nt tracks, one has the charge qj of the j-th track and its rapidity rjj. Given these, one may construct the power 
spectrum of the charge 



P(r/,fc) 



N t 

3=1 



(22) 



Since the qj are conserved numbers and not their densities, one has P(rf, k) = r|P(r/, k). One may similarly construct 
the power spectrum of any conserved quantity: S, B or even the proton number, by replacing qj in the above formula 
by the quantum number under discussion. Again, a plot of P(r/,fc), averaged over events, as a function of k will 
smear out many of the effects that we have discussed. Some information can be gained by comparing the results 
with those obtained from 'mixed events', i.e., when tracks from different events are randomly thrown together into 
an artificial event with the right multiplicity. 

We find that the distinction between Kelly and Fick theories can be observed most directly in the event-to-event 



distribution of r/|n(r/,fc)| oc yP(r/,fc). The proportionality constant requires knowledge of the initial volume. 
However, both the first and second order diffusion equations being linear, the overall normalization of the n and 



v are immaterial. We take advantage of this to work with the dimensionless quantities y P(r/,fc) = rn(r, k)/ao 
and tv(t, k)/ao, where the dimensional quantity ctq need not be specified except when making an actual connection 



to models of the initial state. Even so, a prediction of the distribution of y P(t/, k) is not possible since the initial 
distribution is not known. However, starting from the same initial distribution the two theories give rise to completely 
different distributions at freezeout, and the differences can be investigated. 

An example is shown in Figure [5] Here we started with initial conditions drawn from a Gaussian distribution with 



zero mean and unit variance. Then from eq. (jlip it is clear that for Fick diffusion, the distribution of y P(t/, k) at 

freezeout is also a Gaussian, with standard deviation exp[— Z?/c 2 (1/to — 1/r/)]. Note the strong dependence of the 
variance on k, which is also apparent in Figure [5J The distribution of this quantity in the Kelly theory cannot be 
derived so easily. We determined it through a Monte Carlo simulation, deriving the distribution at freezeout time by 
evolving many different samples of initial conditions. The resulting distribution is clearly different. It does not peak 
at zero, the position of the peak shifts with k, and it has a very long tail. These features are intimately connected with 
the transient amplification phenomenon which has been discussed above. Beyond the range shown in the figures, as 
k decreases, the frequency of small y/P increase, the tail shrinks, and the distinction between Fick and Kelly theories 
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is lost in the limit of k — > 0. Most importantly, the fact that in all theories the long-time behaviour is dominated by 
Bjorken attenuation means that the distribution of P(r/, k) (or its square root) is independent of Tt. Event-to-event 
distributions of this kind therefore seem to be the most promising observable distinction between first and second 
order diffusion. 

One of the most widely studied signals of a critical end-point of QCD is the distribution of event-to-event fluctu- 
ations of conserved quantities contained within a rapidity acceptance window. The present work, with its removable 
simplifications, shows that these measures could be influenced, possibly strongly, by the nature of the transport 
process — whether first or second order. Therefore our understanding of the nature of diffusion needs to be improved 
before event-to-event fluctuations can be interpreted. We have shown that one way to do this is to study the event- 
to-event distribution of the power spectrum of the conserved charge (see eq. [22]) , at tj for different k (see Figure [5]) . 
There are qualitative differences between the results for the two theories of diffusion. 

In summary, here are our main conclusions — 

1. The Fourier coefficients, n{r, k), and the power spectrum, P(r,k), of the number density profile (see eq. [3]), 
obtained by Fourier transforming in the rapidity (which is equal to rj for boost invariant flow) are of interest in 
the study of diffusion and hydrodynamics. One can construct number density profiles for any quantum number 
which is conserved in strong interactions. The obvious ones are B (or N p ), Q, and S. A power spectrum can 
also be constructed from experimental data (see eq. I22p , and easily compared to theory. 

2. The net conserved number is obviously independent of the kind of hydrodynamics. As a result, the event-to-event 
distribution of this quantity tells us directly about initial conditions. The net charge within a given rapidity 
window, however, may evolve with the hydrodynamics (see Figure [4]). Since freezeout may occur at different 
times for different centrality, this may show up as a dependence of the net observed charge on centrality. 

3. The Fourier modes give interesting information on the nature of the transport theory. In Kelly (second order) 
diffusion, the power in short range modes, i.e., modes with large k, could be large. This is impossible in Fick 
(first order) hydrodynamics. The simplest of observables, i.e., a plot of P(t/, k) — t?|ti(t/, k)\ 2 as a function of 
k, for an event, can potentially distinguish between these two theories. One thing to note is that this must be 
done on an event-to-event basis. Averaging over events before constructing the power spectrum could wash out 
the signal. 

4. Event-to-event distributions of the observed power spectrum (see eq. |2"2"1) . or its square root, for one or more 
values of k can potentially distinguish between Fick and Kelly theories of diffusion (see Figure [5] and the 
corresponding discussion in the text). 

5. If the Kelly theory is ruled out by such observables, then one can construct an upper bound on the transport 
coefficient, 

V < r/sinh 2 A?7, (23) 
using the rapidity interval Arj over which all structure has been washed out and an independent estimate of t/. 

This work was motivated by discussions at the INT program "The QCD Critical Point" . SG would like to thank 
the participants, especially V. Koch, G. Roland, M. Stephanov and N. Xu, for discussions, and gratefully acknowledge 
the hospitality of the Institute of Nuclear Theory at the University of Washington during this period. 
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